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A simulated annealing algorithm is employed to generate a stochastic model for a Berea and 
a Fontainebleau sandstone with prescribed two-point probability function, lineal path function, 
and "pore size" distribution function, respectively. We find that the temperature decrease of the 
annealing has to be rather quick to yield isotropic and percolating configurations. A comparison 
of simple morphological quantities indicates good agreement between the reconstructions and the 
original sandstones. Also, the mean survival time of a random walker in the pore space is reproduced 
with good accuracy. However, a more detailed investigation by means of local porosity theory shows 
that there may be significant differences of the geometrical connectivity between the reconstructed 
and the experimental samples. 

PACS: 61.43Gt (Powders, porous materials) 
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I. INTRODUCTION 

The microstructure of porous media determines their 
macroscopic physical properties such as conduc- 

tivity, elastic constants, relaxation times, permeabilities 
or thermal properties. The relation between geometric 
microstructure and physical properties is a fundamental 
open problem whose solution is important to many appli- 
cations ranging from geophysics to polymer physics and 
material science. For geometric modeling of porous me- 
dia it is therefore important to characterize microstruc- 
tures quantitatively and to construct models using given 
geometric characteristics. 

The reconstruction of random media with given 
stochastic properties is also of great interest for a va- 
riety of other reasons, (i) Digitized 3D geometries of 
real sandstones are difficult to obtain. Therefore, the re- 
construction provides a method for easily generating de- 
tailed geometries as needed e.g. in numerical calculations 
of macroscopic material parameters like those mentioned 
above, (ii) The reconstruction of 3D samples from 2D 
data, (iii) Any calculation of macroscopic quantities of 
random media needs a set of stochastic functions which 
describes the geometry. The reconstruction can help to 
decide which functions one should use. The present work 
will focus mainly on the latter question, the so-called in- 
verse problem. 

Recently, a simulated annealing algorithm for the re- 
construction of random porous media with predefined 
stochastic functions was proposed This method was 
used for the reconstruction of two-phase porous media, 
i.e. sandstones, with given two-point probability function 
and lineal-path function, which were measured from real 
sandstones [^,^. In the present article, we will continue 
this work. An advantage of the simulated annealing is 
that it allows the reconstruction of a variety of different 
stochastic functions, with the available CPU time being 
the only limit on the use of a set of functions. In Refs. 

the two-point probability function and the lineal 
path function were used. Here we will extend the in- 
vestigations to the reconstruction of the so called "pore 
size" distribution [0J|], which is not the usual quantity 
obtained from mercury porosimetry. We will study a 
Berea sandstone and a Fontainebleau sandstone. 

A comparison of the reconstructions with the original 
sandstones shows to which extent the characteristics of 
the original geometry are reconstructed. A salient feature 
of the original sandstone is the high degree of connectiv- 
ity of the pore space. Hence, one important criterion for 
judging a stochastic reconstruction method is its ability 
to reproduce the connectivity of the original sandstone. 
It is shown that Berea reconstruction does a better job 
of capturing connectivity than does the Fontainebleau re- 
construction, although neither captures connectivity par- 
ticularly well. Another important criterion is the ability 
of the reconstruction procedure to reproduce the macro- 
scopic properties of the original materials. We show that 



the mean survival times of our reconstructions of both 
sandstones agree well with mean survival times of the 
original sandstones. 

The paper is organized as follows. In section II we 
describe the reconstruction algorithm. In section III we 
introduce the quantities we use in the reconstruction pro- 
cess and for characterizing the sandstones. Because de- 
tailed discussions can be found in the references, we will 
focus mainly on practical aspects and details of the im- 
plementation. In section IV we present the results for re- 
constructions of a Berea and a Fontainebleau sandstone, 
respectively. 

II. THE RECONSTRUCTION METHOD 

A two-phase porous medium consists of a pore or void 
phase P and a matrix or rock phase M. Its microgeometry 
is described in detail by the characteristic function 

r for f e M 
Xp(x) = |^ for feP ■ W 

For a discretized sample x = (jia, ^20, is a) is the position 
vector of a cubic grid with ji — 0, . . . , — 1, lattice 
constant a and size Mi x M2 x M3. The total number of 
grid points is given by = M1M2M3. In the following 
we will refer to a grid point also as a voxel. 

The porosity (/), the probability of finding a point in 
pore space, is given by 

(/.= (Xp(af)). (2) 

With the assumption of a homogeneous, stationary, and 
ergodic stochastic porous medium the angular brackets 
denote a volume average. 

The reconstruction is carried out by means of a sim- 
ulated annealing method Q with target or "energy" - 
function defined as 

X 

The function /rof is the stochastic function to be recon- 
structed, whereas ft is the actual value of this function 
measured at iteration step t. For the reconstruction of 
more than one function the energy Et at the iteration 

step t is given by Et = J2k -^tift'^^) where the index k 
numbers the different functions to be reconstructed. 

Starting from a random configuration with porosity 
(j) two voxels of different phase are exchanged at each 
iteration step. Thus, the porosity (jj remains constant 
during the reconstruction process. The new configuration 
is accepted with the probability given by the Metropolis 
rule 

fl if Et<Et^i 

[e T if > Et-1 
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where T plays the role of a temperature. In the case 
of rejection, the old configuration is restored. By de- 
creasing the temperature T configurations with minimal 
energy E, i.e. with minimal deviations of the stochas- 



tic functions fj;''^ from their reference functions f^gf, are 
generated. The process terminates after a certain num- 
ber of consecutive rejections. Here, the reconstruction 
was finished after 10^ consecutive rejections. 



(fc) 



III. MEASURED QUANTITIES 

A. The two-point probability function 

Using Equation the two-point probability function 
is defined as 



82(^1, ^2) = (Xp(fl)Xp(x2)). 



(5) 



For a homogeneous and isotropic medium S2(25'i,X2) = 
S2(r) with r = \xi — X2\ holds. In this case S2(?') can be 
evaluated without loss of information from the intersec- 
tion of the sample with a plane or even a line. To speed 
up the numerical evaluation of S2 it is therefore sufficient 
to sample S2 only in directions of the principal axis given 
by the unit vectors . 

Hence, S2 is calculated by evaluating equation (j^) 
for every pair of voxels at xi and £2 — xi + rci with 
r = 0, 1, . . . , Tc where is a cut off value determined by 
the system size or a multiple of the correlation length. 
During the reconstruction only those terms in S2 are up- 
dated that have changed due to the exchange of voxels. 
On the sample boundaries we impose periodic boundary 
conditions. 

For reconstruction purposes this simplification, i.e. the 
reconstruction of S2 only in the direction of the three 
coordinate axes, may create some problems. Because 
all other directions remain unoptimized, S2 measured 
in these directions may differ from the reference func- 
tion. If this happens the reconstructed sample is no 
longer isotropic |9|,|l0f|. Furthermore, the reconstruction 
of S2 along orthogonal lines only, reduces the three- 
dimensional optimization problem effectively to the opti- 
mization of three one-dimensional two-point probability 
functions. This may reduce the number of conditions 
that the two-point probability function has to fulfill to 
be realizable and may lead one to conclude incor- 
rectly that the reconstruction is realizable in the space 
dimension of interest. However, the problem of realiz- 
ability does not apply to the reconstructions presented 
here, because the two-point probability functions used 
as reference functions are measured from digitized three- 
dimensional images of real sandstones. Nevertheless, the 
reconstructions have to be checked for their isotropy. 

The isotropy of the reconstructions can be improved 
e.g. by a full reconstruction of the two-point probabil- 
ity function using Fourier transform techniques or by ro- 
tating the sample during the reconstruction. This was 



done for two-dimensional reconstructions in |T^. How- 
ever, for three-dimensional reconstructions this leads to 
prohibitive increase of computation time. 

For the specific surface s, i.e. the surface per unit 
volume of the interface between pore space and matrix 
space, it is known that [l^ 



s = -4 



9S2(r) 



dr 



(6) 



Therefore, a reconstruction of S2 implies that the spe- 
cific surface area of the reconstructions matches that of 
the reference sample. Equation (^) may also be used for 
the calculation of the specific surface. Here, we use a dif- 
ferent, numerically very efficient method introduced in 
|l3| for calculating s. 



B. Lineal path function and "pore size" distribution 

The lineal path function L(r) is defined to be the proba- 
bility of finding a line segment of length r entirely in pore 
space when the line segment is randomly thrown into the 
porous medium Hence, for r = 0, L(0) = (j) holds. 
The lineal path function is related to the linear contact 
distribution introduced in mathematical stochastic ge- 
ometry [^-0 . The lineal path function is calculated as 
follows: For a given pore voxel the lineal path r in units 
of the resolution a is given as the number of pore voxels 
lying between the given pore voxel and the nearest ma- 
trix voxel in direction e^. Evaluating the line segments 
starting from each pore voxel in all three coordinate di- 
rections Ci and counting the number l{r) of line segments 
with length r, the lineal path function is given by 

L(r) = l{r)/ ((Ml - r)M2M^ 

+Mi{M2 - r)K'h + MiM2(M3 - r)) (7) 

where we assume non periodic boundary conditions. 

The lineal path function incorporates information 
about the connectivity of the pore space. Of course it 
would be desirable to reconstruct functions which pro- 
vide a more complete information about the geometric 
connectivity of the pore space, as for example the cluster 
correlation function |^] or the local percolation proba- 
bility 1^, but currently the computation time for eval- 
uating these functions prevents their use in the above 
reconstruction scheme. 

The "pore size" distribution function P{5) is defined 
such that P((5)d(5 is the probability that a randomly cho- 
sen point in the pore space lies at a distance [(5, 5 -\- dS] 
from the nearest point on the interface 
to the spherical contact distribution |1 
sociated cumulative distribution function 




It is related 
111. The as- 



P{5)d5 



(8) 
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gives the fraction of the pore space which has a diameter 
larger then So- Clearly 



and 



P(0) = 



df 

36 



(9) 



(10) 



The mean pore size is given as 

{S)= S P{d)dS ^ / F(5)d(5. (11) 
Jo Jq 

The quantity P(S) arises in rigorous bounds on the mean 
survival time [p[. 

We compute an approximation to P{5) by choosing a 
random point in pore space and measuring its distance 
i5 to the nearest point on the matrix-pore interface as- 
suming periodic boundary conditions. This process is 
repeated for several random points in pore space. The 
"pore size" distribution is then obtained by binning the 
distances S and dividing by the number of random place- 
ments in pore space. We emphasize that the random 
placements in pore space are not necessarily grid points. 
The computation of 6 is only approximate because it re- 
quires a modeling of the interface between pore and ma- 
trix space. Here, we assume that the internal surface 
is given by the surface of the cubic voxels. This is the 
same modeling which is used e.g. in a computation of 
the mean survival time or in finite difference calculations 
of transport properties. In general this may overestimate 
the specific surface area appearing in Equation ( p^ by a 
factor of roughly 1.5 | p^ . 

For use in the reconstruction we measured 6 as the 
distance between a pore voxel and the nearest matrix 
voxel. The resulting function is not equal to P{S) because 
i5 can now only take values of S — y^i^ + p + A;^ with 
i,j,k e Z. 



C. The total fraction of percolating cells 

The total fraction of percolating cells p is a key mea- 
sure in local porosity theory [^,^. Local porosity the- 
ory measures the fluctuations of morphological quanti- 
ties, e.g. porosity, specific surface, connectivity, in cubic 
subsamples of the total sample |2^. We will refer to 
such a cubic subsample with side length L as a mea- 
surement cell. Based on the scale-dependent morpholog- 
ical quantities, the theory provides scale-dependent esti- 
mates for transport parameters from a generalized effec- 
tive medium theory [^ , ^ . 

Let /i(0, L)d(/) be the probability that a given mea- 
surement cell of side length L has a porosity in the in- 
terval [0, (/) -|- d(j)\. The probability density function is 
called local porosity distribution. The probability that 



a given measurement cell with porosity 4> is percolating 
in all three directions is the local percolation probability 
function X{(j),L). Here, percolating in all three direc- 
tions means that each face of the measurement cell is 
connected to the opposite face with a path lying entirely 
in pore space. Using this, the total fraction of percolating 
cells is given by 



p{L)= / fi{cb,L)\{cl>,L)dcb. 



(12) 



The total fraction of percolating cells p{L) is the proba- 
bility of finding a measurement cell with side length L, 
which is percolating in all three directions. Hence, p is a 
measure for the geometrical connectivity. 



D. The mean survival time 

Contrary to the previous quantities the mean survival 
time r is not a purely geometrical but a physical ob- 
servable. The mean survival time r is the average life 
time of a random walker, which can freely move in pore 
space, but gets instantly absorbed on contact with the 
pore-matrix interface. It is a measure of a characteristic 
pore size. The mean survival time r is calculated using 
a first passage cube (FPC) algorithm |l|,|2|. The FPC 
algorithm uses the fact that the mean time it takes for 
a certain type of diffusive random walker starting at the 
center of a cube with side length 2L to cross the surface 
of this cube is given by 



t{L) « 0.22485X2 



(13) 



Hence, it is not necessary to simulate the steps of the 
walker in detail. Instead one determines the biggest cube 
centered around the position of the walker, which is still 
entirely in pore space. The walker then jumps on the 
surface of this FPC and a time given by Equation ( [T^ ) 
is added to its life time. This procedure is iterated un- 
til the walker touches the interface and gets absorbed. 
The mean survival time is given by averaging over many 
walkers. 

The probability with which the walker jumps to a cer- 
tain point on the surface of the FPC is described by a 
probability density function w{y, z) where y and z are 
the coordinates on the surface assuming without loss of 
generality that x = ±L. For an analytic expression of 
w{y, z) we refer to 



IV. RESULTS 

In this section we present results for reconstructions 
of a Berea sandstone and a Fontainebleau sandstone. 
For both sandstones we computed reconstructions with 
the lineal path function and the two-point probabil- 
ity function (LS reconstruction) as well as reconstruc- 
tions with the "pore size" distribution function and the 



4 



two-point probability function (PS reconstruction). The 
Berea sandstone and its reconstructions have dimension 
128 X 128 X 128 and resolution a = lO/xm. The porosity 
is (j) = 0.1775. The Fontainebleau sandstone has dimen- 
sions 299 X 300 X 300 its reconstructions have dimension 
128 X 128 X 128. The resolution is a = 7.5//m, the porosity 
is (f> — 0.1355. The reconstructed functions were calcu- 
lated as described above, i.e. we used periodic boundary 
conditions except for the lineal path function. The two- 
point probability function S2 (r) was reconstructed in the 
interval r — 0, 1, ... , 63. The annealing process termi- 
nated after 10^ subsequent rejections. 

We performed 5 LS reconstructions and 5 PS recon- 
structions for both the Berea and the Fontainebleau 
sandstone. Some of the results are summarized in Ta- 
ble I for the Berea and Table || for the Fontainebleau 
sandstone, respectively. The values are averaged over 
5 reconstructions. The quantity tD is the mean sur- 
vival time multiplied by the diffusion coefficient D for 
the random-walk process discussed above. The quantity 
fp gives the fraction of pore voxels which belong to the 
percolating cluster. 

Figure |] shows 2D slices of the original and the re- 
constructed sandstones. The top row shows the origi- 
nal sandstones, the row in the middle the LS reconstruc- 
tions, and the bottom row the PS reconstructions. The 
slices were taken from reconstructions with the value of 
p{L = 60) close to the average values given in the tables. 
All slices are chosen to have average porosity. In the 
case of the Berea sandstone the reconstructions look sim- 
ilar to the original sandstone while for the Fontainebleau 
sandstone the reconstructions are clearly distinguishable 
from the original sandstone. The matrix of the origi- 
nal Fontainebleau sandstone shows a granular structure 
where single grains can be identified. The pores between 
these grains are long and narrow. In the reconstructions 
no granular structure of the matrix space is visible. The 
pores of the reconstructions are more rounded in shape. 
For both sandstones the number of isolated pores is sig- 
nificantly higher in the reconstructions. This is expressed 
by fp, the fraction of pore space belonging to the perco- 
lating cluster. For the Berea 97.16% of the pore space 
belongs to the percolating cluster whereas for the recon- 
structions this fraction is roughly 10% smaller. For the 
Fontainebleau the difference is even bigger. In the origi- 
nal sandstone 99.35% of the pore space percolates while 
for the reconstructions fp is approximately 52%. Here 
we find that one LS and one PS reconstruction of the 
Fontainebleau is not percolating in all three directions. 

In the course of our work we used also a slowly de- 
creasing step function for the temperature T to obtain 
an optimal match of the reconstructed functions. Sur- 
prisingly, with a slow cooling schedule the majority of 
the reconstructed configurations was not percolating in 
all three directions. Furthermore, the reconstructed sam- 
ples showed a strong anisotropy with £{$2) of order 10~^ 
where S2 denotes the two-point probability function mea- 



sured in the directions -I- e^^i . 

For the reconstructions presented here, we used a fast 
exponential cooling schedule T — exp (j^) where t as 
above denotes the iteration step. This cooling schedule 
took approximately 30 • N iteration steps to complete a 
reconstruction, whereas the slow cooling took more than 
300 • N iterations steps. Using the fast cooling schedule 
all reconstructions of the Berea sandstone are percolating 
in all three directions and only one LS and one PS recon- 
struction of the Fontainebleau is not percolating. Also 
with the fast cooling schedule the reconstructed func- 
tions are matched very well, i.e. E{S2) is of order 10"^", 
E{L) is of order 10~*. Moreover, the anisotropy mea- 
sured in terms of E{S2) was reduces by an order of mag- 
nitude. Plotting S2 only the PS reconstructions of the 
Fontainebleau showed small deviations. 

Our explanation for the fact that a slower cool- 
ing schedule results in reconstructions with stronger 
anisotropy and only poor connectivity is the artificial 
anisotropy introduced by reconstructing S2 and L only 
in three directions. With increasing number of iterations 
the influence of the isotropic, random starting configura- 
tion is decreased while the anisotropic calculation scheme 
of the two-point probability functions as described above 
becomes more significant. This view agrees with pre- 
vious work where 3D isotropic reconstructions of 
sandstones using also a fast cooling schedule were pre- 
sented. The poor connectivity of the reconstructions 
with slow cooling schedule may be a result of their strong 
anisotropy. 

In Figures ^, ^ and ^ the two-point probability func- 
tions S2, the lineal path functions L, and the "pore size" 
distribution function P, respectively, are plotted for both 
sandstones using lines and for typical reconstructions us- 
ing dots. Here, typical means that the energy given by 
Equation (^) of the reconstructed functions is close to 
the average value. In the case of the two-point proba- 
bility function, the reconstructed functions appear to be 
indistinguishable from the reference functions. The same 
applies to the lineal path functions measured from the 
LS reconstructions. The lineal path functions of the PS 
reconstructions clearly underestimate the reference func- 
tions. Complementary to this, P is equally matched by 
both types of reconstructions as shown in Figure ^. A 
logarithmic plot of P and L reveals that the reconstruc- 
tions only poorly match the tails of those functions. This 
is due to the extremely small values of P{6) and L(r) for 
large d and large r, respectively. Furthermore, for the 
lineal path function L this may be a finite size effect be- 
cause L is a long ranged function with L(r) > for values 
of r in the order of a third of the system size. 

Looking at Figure ^ it seems that for our reconstruc- 
tions the lineal path function L and the two-point prob- 
ability function S2 incorporate nearly the same infor- 
mation about the shape of the pores as S2 and P do. 
Moreover, as seen from Figure ^ the PS reconstruction 
lacks information about long line segments. This may be 
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understood from the fact that P is a very short ranged 
function with P{S) 7^ only in a range smaller than the 
correlation length. For the Berea sandstone P{S) — Q for 
6 > 60/im and for the Fontainebleau sandstone P((5)=0 
for S > 78//m. Hence, even though P contains full three- 
dimensional information about spherical regions in pore 
space, similar information is already provided by the two- 
point probability function S2. 

Looking at the parameters given in Table | and g 
the PS reconstructions match better the specific surface. 
This is expected because the specific surface s can be 
measured from either S2 or P as seen from Equation (j^) 
and However, the value of s measured from Equa- 

tion ( p^ turns out to be roughly a factor 1.5 bigger than 
the value computed from Equation (|^). This is due to 
the simple surface modeling in the calculation of P as 
discussed above. We believe the kink in P{6) for 5 — 0.5a 
to be an artifact of the discretization. 

In the case of Berea we find the best agreement of the 
mean pore size (6) for the PS reconstructions, while in 
the case of the Fontainebleau the match appears to be 
better for the LS reconstruction. The lineal path func- 
tion seems to be better suited to describe the long nar- 
row pores of the Fontainebleau sandstone than the "pore 
size" distribution. As already seen from the 2D slices 
the appearance of the Berea sandstone is quite similar 
to the appearance of the reconstructions. The pores are 
much more rounded in shape. Here, the reconstruction 
is slightly improved with respect to the mean pore size 
by incorporating P, which contains information about 
spherical regions. 

We find analogous results for the mean survival time 
T, which is a diffusive transport property. In fact the 
mean survival time can be related to the mean "pore 
size" [|j . The mean survival time is a physical transport 
property, but unlike the fluid permeability, it does not 
capture information about the dynamical connectivity of 
the pore space; indeed, neither does the conductivity (or 
formation factor) . Nonetheless, a cross-property formula 
relating the fluid permeability to a combination of the 
porosity, mean survival time and formation factor F has 
been shown to be a highly accurate estimate of k for 
sandstones [^. This cross-property formula was used 
to demonstrate that the permeability of another recon- 
structed sandstone |^ was in excellent agreement with 
the exact Stokes solution determination of the perme- 
baility of the original sandstone . 

We note that a good match of r alone may not al- 
ways indicate a good match for J- or k. In fact, com- 
bining results from |^,^ and |^ suggest that T cor- 
relates strongly with the local percolation probability p 
which is a measure of the geometrical connectivity. We 
find significant differences between the local percolation 
probabilities of the real sandstones and the present re- 
constructions. 

Figure || shows plots of the local percolation probabil- 
ity p for the original sandstones, the LS reconstructions, 
the PS reconstructions, and reconstructions of the two- 



point probability function only (S reconstruction). The 
pure S reconstructions are included here for comparison 
to previous work The curves shown for the recon- 
structions are averaged over five configurations each. The 
local percolation probability p of the reconstructions lie 
well below the curves of the original sandstones. From 
this plot it seems, that neither the use of L nor the use 
of P can significantly improve the geometric connectiv- 
ity compared to the reconstruction of S2 only. The dif- 
ferences between the curves of the three reconstructions 
seems to be within the range of statistical fluctuations. 
Nevertheless, a similar result for the S reconstruction of 
the Fontainebleau sandstone was presented in for a 
larger sample. Other work j|] showed that a LS recon- 
struction of a different Fontainebleau sandstone repro- 
duced the geometric connectivity well. 

In general comparing the LS reconstructions and the 
PS reconstructions the resulting configurations are quite 
similar. For our reconstructions the two-point proba- 
bility function S2 and the lineal path function L incor- 
porate nearly the same morphological information as S2 
and P do. Looking at the reconstructions of the Berea 
sandstone and the reconstructions of the Fontainebleau 
sandstone it appears that the latter one is much more 
demanding to reconstruct. This may be due to its char- 
acteristic granular structure, the narrow pore throats, the 
lower porosity, and to the larger sample size of the orig- 
inal sandstone. We also note that the Berea sample is 
only 128 x 128 x 128 resulting in poor statistical quality. 

Our work has shown that simulated annealing provides 
a flexible and simple to implement method for recon- 
structing two-phase random media and that local poros- 
ity theory provides highly sensitive tools for their com- 
parison and analysis. However, with present computer 
power there is still a need to introduce simpliflcations to 
reduce the computation time. Reconstructing the two- 
point probability function S2 only in certain directions, 
may introduce artiflcially a strong anisotropy or affect 
the connectivity. We find, that a fast cooling schedule 
can reduce this problem. This implies that the final con- 
figuration is not completely independent from the initial 
configuration, and hence the reconstructed microstruc- 
ture does not only depend on the reconstructed statistical 
functions as would be desirable. 
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TABLES 





Berea 


LS 


PS 


porosity 


0.1775 


0.1775 


0.1775 


specific surface [mm~^] 


13.9 


14.7 


14.4 


mean survival time tD [^m^] 


100 


89 


93 


mean pore size {5} [/jm] 


6.71 


6.52 


6.66 


/p [%] 


97.16 


88.76 


85.49 


p{L = 60) 


0.997 


0.747 


0.712 



TABLE I. Characteristic quantities of the Berea sandstone, its LS reconstructions, and its PS reconstructions. The values 
given for the reconstructions are averaged over five configurations. 





Fntbl 


LS 


PS 


porosity 


0.1355 


0.1355 


0.1355 


specific surface [mm^^] 


10.0 


10.6 


10.4 


mean survival time tD [fim^] 


134 


129 


121 


mean pore size {6} [/um] 


7.85 


7.88 


7.73 


fp [%] 


99.35 


52.22 


51.26 


p(L = GO) 


o.ono 


0.265 


0.234 



TABLE II. Characteristic quantities of the Fontainebleau sandstone, its LS reconstructions, and its PS reconstructions. The 
values given for the reconstructions are averaged over five configurations. 
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FIGURES 



see figure_l.jpg 

FIG. 1. Two dimensional slices of the Berea sandstone, the Fontainebleau sandstone, and two reconstructions of each. The 
top row shows the original sandstones, the middle row the LS reconstructions, and the bottom row the PS reconstructions with 
the Berea sandstones on the left and the Fontainebleau sandstones on the right side. All slices have approximately the average 
porosity = 0.1775 for the Berea and 4> = 0.1355 for the Fontainebleau, respectively. 



Berea 
Berea LS 
Berea PS 
Fntbl 
Fntbl LS 
Fntbl PS 




300 



FIG. 2. Two-point probability functions S2 of the Berea (top) and the Fontainebleau (bottom) sandstone. The solid lines 
show the reference functions, the points show typical reconstructed functions for a LS and a PS reconstruction. 



9 




10 




11 




10 20 30 40 50 60 

L [in units of a] 




L [in units of a] 

FIG. 5. Total fraction of percolating cells p for the Berea sandstone (top) and the Fontainebleau sandstone (bottom). The 
solid lines show p measured from the original sandstones. The curves shown for the reconstructions are averaged over five 
configurations each. 
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